Exactly solvable models through the generalized empty in- 
terval method, for multi-species interactions 



A. Aghamohammadi'^''^^, M. Alimohammadi^^, & M. Khorrami'^'^'^ 

Department of Physics, Alzahra University, Tehran 19834, Iran. 
^ Physics Department, University of Tehran, North Karegar Avenue, Tehran, 

Iran. 

Institute for Advanced Studies in Basic Sciences, P. O. Box 159, 

Zanjan 45195, Iran. 
Institute of Applied Physics, P. 0. Box 5878 Tehran 15875, Iran. 

Abstract 

Multi-species reaction-diffusion systems, with nearest-neighbor interac- 
tion on a one-dimensional lattice are considered. Necessary and suffi- 
cient constraints on the interaction rates are obtained, that guarantee 
the closedness of the time evolution equation for E^{t)'s, the expectation 
value of the product of certain linear combination of the number oper- 
ators on n consecutive sites at time t. The constraints are solved for 
the single-species left-right-symmetric systems. Also, examples of multi- 
species system for which the evolution equations of E'^(t)'s are closed, are 
given. 

PACS numbers: 05.40.-a, 02.50.Ga 
Keywords: reaction-diffusion, generalized empty-interval method, 

multi-species 



^ mohamadi@azzahra. ac. ir 

^alimohmdCSut.ac.ir 

•^mamwadCiiasbs.ac.ir 



1 Introduction 



There is a well-established framework for equilibrium statistical mechanics, but 
thermal equilibrium is a special case, and there isn't a corresponding straight- 
forward framework for investigating the properties of systems not in equilib- 
rium. There is no general approach to systems far from equilibrium. Different 
methods have been used to study these models. These include analytical and 
asymptotic methods, mean-field methods, and large-scale numerical methods. 
For high-dimensional systems, mean-field techniques give exact or reasonable 
approximate results. But their results for low-dimensional systems are gener- 
ally not adequate. So, people are motivated to study stochastic models in low 
dimensions. Models in low dimensions, should also be in principle easier to 
study. Exact results for some models on a one-dimensional lattice have been 
obtained, for example in .1..2...3..4 5. 6, 7, 8, 9, 10, 11 , 12. 13, 14, 15, 16..17). 

The term exactly-solvable have been used with different meanings. In |18| . 
a ten-parameter family of reaction-diffusion processes was introduced for which 
the evolution equation of n-point functions contains only n- or less- point func- 
tions. The average particle-number in each site has been obtained exactly 
for these models. In T^'TU^, the same method has been used to analyze the 
above mentioned ten-parameter family model on a finite lattice with boundaries, 
and in |21j similar method has been used to study models with next-nearest- 
neighbor-interactions. In 122] and integrability means that the A^-particle 
conditional probabilities' S-matrix is factorized into a product of 2-particle S- 
matrices. Another method which has been used to solve some reaction diffusion 
models exactly is the empty interval method, and its generalizations. 

The empty interval method (EIM) has been used to analyze the one di- 
mensional dynamics of diffusion- limited coalescence |231|2IESlEni • Using this 
method, the probability that n consecutive sites are empty has been calculated. 
For the cases of finite reaction-rates, some approximate solutions have been 
obtained. EIM has been also generalized to study the kinetics of the g-state 
one-dimensional Potts model in the zero-temperature limit |27j . 

In j28) , all the one dimensional reaction-diffusion models with nearest neigh- 
bor interactions which can be exactly solved by EIM have been studied. EIM 
has also been used to study a model with next nearest neighbor interaction |29) . 
In |30|. exactly solvable models through the empty- interval method, for more- 
than- two-site interactions were studied. In the conventional EIM has been 
extended to a more generalized form. Using this extended version, a model 
which can not solved by conventional EIM has been studied. 

In this article, we consider systems, in them particles of more than one 
species interact on a one-dimensional lattice. The interaction is nearest-neighbor. 
Each site of the lattice, either is empty, or is occupies by one particle. In section 
2, we seek necessary and sufficient conditions on the reaction rates, so that the 
time evolution equation for E^^{t) is closed. This quantity is the expectation 
of the product of a specific linear combination of the number operators (cor- 
responding to different species) at n consecutive sites beginning from the fc-th 
site. In section 3, all single-species left-right symmetric reaction-diffusion sys- 
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terns solvable through the generalized empty-interval method, are classified. In 
section 4, multi-species systems are investigated, which are solvable through the 
generalized empty-interval method, but are effectively single-species. In section 
5, some specific families of two-species systems are investigated, which are ex- 
actly solvable through the generalized empty-interval method. Finally, section 
6 is devoted to concluding remarks. 

2 Models solvable through the generalized empty- 
interval method 

To fix notations, let us briefly introduce the multi-species reaction-diffusion 
systems with nearest-neighbor interactions, on a periodic lattice. Let the lattice 
have L + 1 sites. The observables of such a system are the operators N", where 
i with 1 < i < L + 1 denotes the site number, and a with 1 < a < p -I- 1 denotes 
the type of the particle. One can regard a = p+1 as a vacancy. A^f is equal to 
one, if the site i is occupied by a particle of type a. Otherwise, N" is zero. We 
also have a constraint 

SaiVf = 1, (1) 

where s is a covector the components of which (sq's) are all equal to one. The 
constraint simply says that every site, either is occupied by a particle of 
one type, or is empty. A representation for these observables is 

AT^ 1 ® . . . ® 1 07V" 1 ® • • • (g) 1, (2) 

i-l L+l~i 

where N" is a diagonal (p -I- 1) x (p+1) matrix the only nonzero element of 
which is the a'th diagonal element, and the operators 1 in the above expression 
are also (p+1) x (p+1) matrices. It is seen that the constraint ^ can be 
written as 

s • N = 1, (3) 

where N is a vector the components of which are N°"s. The state of the system 
is characterized by a vector 

P G V ® • • • (g) V, (4) 

L + l 

where V is a (p + l)-dimensional vector space. All the elements of the vector P 
are nonnegative, and 

S ■ P = 1. (5) 

Here S is the tensor-product of i + 1 covectors s. 

As the number operators N" are zero or one (and hence idempotent), the 
most general observable of such a system is the product of some of these number 
operators, or a sum of such terms. 
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The evolution of the state of the system is given by 



P = 7^ P, (6) 

where the Hamihonian Ti is stochastic, by which it is meant that its nondiagonal 
elements are nonnegative and 

S 7^ = 0. (7) 
The interaction is nearest-neighbor, if the Hamiltonian is of the form 

L+l 

= ^ Hi,i+i, (8) 

where 

:= l(8)--(gil (8ig(g) l(g)- - -(8)l . (9) 

i—l L—i 

(It has been assumed that the sites of the system are identical, that is, the 
system is translation-invariant. Otherwise H in the right-hand side of (jS)) would 
depend on i.) The two-site Hamiltonian H is stochastic, that is, its non-diagonal 
elements are nonnegative, and the sum of the elements of each of its columns 
vanishes: 

(s®s)i7 = 0. (10) 
Now consider a certain class of such observables, namely 

fe+n-l 

5,-,„:= n (a. NO, (11) 

where a is a specific {p+ l)-dimcnsional covector, and N.^ is a vector the compo- 
nents of which are the operators N". We want to find criteria for H, so that the 
evolutions of the expectations of £^ „ 's are closed, that is the time-derivative of 
their expectation is expressible in terms of the expectations of £'^„'s themselves. 
Denoting the expectations of these observables by „ , 



El„ := S ffc«„P, (12) 



we have 



+ ^ S £^^^Hk-i+i,k+i P 
1=1 

+ S £^ nHk+n-l,k+n P- (13) 



3 



From this, and using ||2Jl, it is seen that the necessary and sufficient conditions 
that the left-hand side be expressible in terms of „ 's are 

(s ® s)[(s • N) (a • N)]i? =^l(s ® s)[(s • N) ® (a • N)] 

+ 0l(s®s)[(s-N)® (s-N)] 
+ i/L(s®s)[(a-N)(g)(a-N)], (14) 
(s s)[(a • N) (a • N)]H =A(s s)[(a • N) ® (a • N)], (15) 
(s ® s)[(a • N) ® (s • N)]H =^r(s ® s)[(a • N) ® (s • N)] 

+ 0R,(s®s)[(s-N)® (s-N)] 

+ ^'R(s®s)[(a-N)® (a-N)], (16) 

for some arbitrary numbers A, /iL,i?,, Ol,r, and Vh^R- Using the identity 

s(b • N) = b, (17) 

for an arbitrary covector b, one arrives at 

(s (g) a)H =/iL(s a) + 9i^{s ® s) + j/L(a ® a), 
(a(g)a)ff =A(a(g)a), 

(a ig) s)i? =/iR(a (g) s) + 6'r(s (g s) + i/R(a (g a). (18) 

The Hamiltonian H should of course satisfy pOfl as well, and its nondiagonal 
elements should be nonnegative. In this case, the evolution equation for E^^, 
for < rt < L + 1, becomes 

^k.n = ^L-Bfc_^i „_l+i^L-E'^_l_„_,.l+6'R£;^„_l+^'R£;^„+l + [/iL+(?^-l)A+/iR]£^fc „. 

(19) 

For n — L + 1, one has 

KL+i = (i + i)A£;iVi- (20) 

For n = 0, from (jSJ we have the boundary condition 

Elo = 1- (21) 

Equations (|19|l and H2U|I . and the boundary condition l|21() are a closed set 
of evolution equations for Efj's. These equations are quite similar to those 
obtained in [2HI- In fact, the case there is a special case of what considered 
here, with p = 1 and a = (0, 1). Although the criterion for the closedness of the 
evolution equations does depend on p and a, the evolution equations for Efj 's do 
not depend on these. In [28', situations were considered in which A was zero, so 
that the evolution of the block comes solely from its ends. This makes solving 
the evolution equations easier. Finally, since the system under consideration 
is translationally-symmetric, if the initial condition is translationally- invariant, 
the state of the system would remain translationally-invariant at all times. In 
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this case, „ does not depend on k, and one would have 



K =(^L + dR)E^_^ + {i^L + i^r)E^+i + [ml + {n~ 1)A + mr)£;^, < n < L + 1, 

=1- (22) 

The system is left-right symmetric, iff the Hamiltonian is invariant under 
permutation: 

nHn = H, (23) 

where 11 is the permutation matrix: 

n(u (g) v) = V (g) u. (24) 
It is easily seen that if H satisfies (|23|) and p8|) . then 

Ml = MR = f R = O-l^Ok^ 0. (25) 

3 Classification of the single-species left-right- 
symmetric reaction-diffusion systems, which 
are solvable through the generalized empty- 
interval method 

For a single-species system, the vector space V is two-dimensional. Take a 
covector a, which is not a multiple of s. (The case a proportional to s is 
trivial, as s • N = 1). The set i? :— {a g) a, a g) s, s g) a, s (g s}, is a basis 
for V (g) V. If _ff satisfies (|18|l (and of course (|10|l ). then one has the matrix 
elements of the Hamiltonian in this basis. However, not every Hamiltonian 
in this form represent a stochastic system. The nondiagonal elements of the 
Hamiltonian in the physical (ordinary) basis should be nonnegative. So, for the 
single-species systems, the task of classifying the systems solvable through the 
generalized empty-interval method, reduces to finding the matrix elements of 
the Hamiltonian in the physical basis (in terms of the covector a, and the scalars 
A, Mi J': E^nd 9); and imposing the criterion that the nondiagonal elements of H 
in the physical basis be nonnegative. Taking the covector a like 

a=(ai,a2), (26) 

and noting that oi ^ 02 (otherwise a would be proportional to s), and that one 
can rescale a without changing the conditions l|18(l , it is seen that one can take 

ai = e + 1, a2 = C - 1, (27) 

without loss of generality. Then, imposing the criterion that the nondiagonal 
elements of the Hamiltonian are nonnegative, and assuming left-right symmetry, 
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one arrives at the foUowing set of inequahties. 

(1 + e')A + 2^2/1 + 2^0 > |2^(A + /i)|, 
-(1 + e')(2Ai + A) - Aiu - 2ie > |2[(1 + + ^(A + 2^i) + e]\, 
-2/1 > (1 - e^)A - 2fn - 2£,0 > \2[{f - + Cm + (28) 

where 

A -A + 2^iy. (29) 

We also define 

A± :=A ± 2z/. (30) 

We have to solve the inequalities H28|) . for a given value of ^. It is seen that 
changing the signs of ^, v, and 9 simultaneously, while keeping the signs of A, 
yU, and A fixed, does not change the inequalities. So it is sufficient to solve the 
inequalities for nonnegative ^. 

The detailed calculations can be found in the appendix. The results are the 
following. 

i) 

^ = 1, /i 6* = 0, iy<0, 4v< X<2iy. (31) 

The reactions for systems in this class are 

^ 00, with rate X - if 

AA 9A, Ad) with rate 2iy - A, (32) 

for which the evolution equation of (ni ■ ■ ■ Uk) is closed, where rii is the number 
operator at the site i. 

ii) 

< ^ < 1, /i<0(/t=-l) (A~,A+) inside the tetragon ABCi?, t satisfies (|^ . 

(33) 

The coordinates of the vertices of the tetragon ABCE are 
M-l^.A B(uP^). cfi^.A. E( ' + l + 



where the first coordinate is A~, and the second is A+. As an example in this 
class, take ^ = 1/2, and A+ — A~ = 1, which lead to a system with following 
interactions. 

« . . , 3 + 461 

A9 ^ 0A, with rate 

AA, 00 0A, A0 with rate 

AA, ^0, 9A 00, with rate 





16 


3 


+ 461 




16 


9 


- 129 




16 


1 


+ 46* 



A^,(l)A-^AA, with rate . (34) 
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For this system, the evolution equation of {[2ni — (1/2)] • • • [2nk — (1/2)]) is 
closed. 



iii) 

1<^, ^^<0{^^ = -l), r = -^±^, A± = (35) 

It is seen that here there is no allowed region for the rates, but a single point 
(apart from scaling) for each value of ^. Systems in this class, correspond to the 
reactions 

Ad), -> 00, with rate ^ + ^ 
A0, 0A ^ AA, with rate ^ ~ (^6) 
and for them the evolution equation of ((2ni + ^ — 1) ■ • • (2nfe + ^ — 1)) is closed, 
iv) 

^ = 0, ^ (^ = —1), (A~, A+) is inside the square AqBoC'o-Do- 

T satisfies (jS7||. (37) 

The vertices of the tetragon AqBqCqDo are 

Ao(-l,l), Bo(l,3), Co(3,l), Do{l,-l). 



As an example in this class, take A+ = A =1, which lead to a system with 
following interactions. 

A0 ^ 0A, with rate - 
4 

AA, 00 0A, A0 with rate i 

^ 2^ 

AA, A0, 0A ^ 00, with rate — - — 

1 -I- 9/9 

00, AID, IDA ^ AA, with rate — (38) 
For this system, the evolution equation of {{2ni — 1) • • • (2nfc — 1)) is closed, 
v) 

^=1, fi^O{fi = -l), (A",A+) is inside §. r satisfies GO)). (39) 

The region §, is the region limited by the lines AiB' , the horizontal line passing 
through B' , and the line passing through Ai with the slope —1, containing the 
point (1,1), where 

Ai(0,l), S'(0,2). 
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As an example in this class, take A+ = A =1, which lead to a system with 
following interactions. 

Q 

A% ^ 0A, with rate - 

Q 

AA, 00 0A, A9 with rate - 

AA, AiJ), 9A 00, with rate 1 - 6. (40) 
For this system, the evolution equation of (ni • • • Uk) is closed. 



4 Effectively single-species systems 

For a specific {p+ l)-dimensional covector, we seek Hamiltonians satisfying p8|l 
with some values of A, h^^r, v-l^r, and 0l,_r- There are cases, however, where 
p-species systems are effectively single-species. Suppose we can decompose the 
states of the system into two subsets 1 and 2. Corresponding to these states, one 
defines the covectors E^, and E'^. E^, for example, is the sum of the covectors 
corresponding to the states belonging to the first subset (the microstates of the 
first state). It is clear that 

E^ + E^ = s. (41) 

To have an effectively two-state system, the probability that the system goes 
from one of the microstates of the state i to the state j, should not depend on 
the microstates. In terms of the Hamiltonian, this means 

W H ^^WjW . (42) 

3 

In this case, the system described by the Hamiltonian iJ is a two-state sys- 
tem. If, moreover, one seeks informations about the original system, which are 
expressible in terms of only the states 1 and 2 (and not dependent on the mi- 
crostates) then the system is an effectively two-state system. For example, if 
the states of a system are white, blue, and red, one can define the states white 
and colored, provided the probability that the system changes from red to white 
is the same as that of changing from blue to white. If in addition, we are only 
interested in probabilities of the system being white or colored, then the system 
is effectively a two-state system. Obvious generalizations of these systems with 
nearest-neighbor interactions on a lattice, are systems for which 

(E' ® E^' ) iJ = ^ if ''^.i (E'^ E' ) . (43) 

k,l 

An example is a system consisting of two kinds of particles {A and B) dif- 
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fusing on a lattice: 

AB ^ BA, with rate A' 

0A, with rate A 

50 ^ 0B, with rate A. (44) 

This system is effectively one-species, as far as one is concerned only with prob- 
abilities of finding particles (not of a specific type). 

Now let us return to the problem of finding solutions to (|18|l , with a specified 
covector a. We want to show that if the components of the covector a take 
only two values, then the system under consideration is an effectively two-state 
system (or an effectively single-species system, with the states occupied and 
empty). If the components of a take only two values, then one can write a as 

a=aiEi+a2E^ (45) 

where E"s are covectors with the property that their components are either 
zero or one, and their sum is equal to s. It is seen that the covectors s and a 
are linear combinations of and E^, and vice versa. So a result of equations 
(O and lUnjl is that (E' ® W)H is a linear combination of (E'' (g) E')'s; that is, 
H43|l holds. Also, one notes that 

a • N = aiE^ • N + azE^ • N, (46) 

which means that the information we seek involves the probabilities correspond- 
ing to only the subsets 1 and 2. So, multi-species systems solvable through the 
generalized empty-interval method, with covectors the components of which 
take only two values, are effectively single-species. 



5 Some two-species examples 

The states of a two-species system on each site of the lattice can be represented 
by A, B, and 0, the latter being a vacancy. As the first example, consider a 
system for which the Hamiltonian is symmetric, which means the rate of each 
reaction is equal to the rate of its reverse reaction. Also let A = in (|18|) . For 
the three-dimensional covector a, take the choice 

a= (1,-1,0). (47) 

One can then solve 118|l and (|10|l for H and /Ul,_r, t^L.Ri and 6'L,i?,. Also the rates 
(the nondiagonal elements of H) should be nonnegative. These constraints lead 
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to a system with the following reaction rates. 



AA ^ BB; 


with rate 2u + v - 


h 2w 


AB ^ BA; 


with rate 27i + u - 


h 2w 


A% ^ 50, 


with rate 2w 




A% ^ 0A, 


with rate 2u 




A% ^ 05, 


with rate 2u 




A0 00, 


with rate 2v 




50 ^ 0A, 


with rate 2u 




50 ^ 05, 


with rate 2u 




50 ^ 00, 


with rate 2v 




%A ^ 05, 


with rate 2g 




0A ^ 00, 


with rate 2r 




05 ^ 00, 


with rate 2r, 





with the condition 



2q = v + 2w. 



For these rates, 



ML = MR = — (4u + 2z; + 4zi;), 
i^L = i^R = = 6'r = 0. 



(48) 
(49) 

(50) 



In this example, the evolution of 5J^'s are very simple. In fact the evolution 
equations decouple and one has 



En - 
El+i = 0, 
= 1. 



< n < L + 1, 



(51) 



The second example is less trivial. Let H satisfy (|23|l . that is the system has 
left-right symmetry. Also let the covector a be 



a= (1,C,0), 



(52) 



and A = 0. From the first (or third) equations of (|18|1 . one can find /i, v, and 9 
in terms of the rates and the parameter ^. The remaining equations are linear 
equations for rates. However, we have inequalities as well, namely the rates 
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should be nonnegative. A specific class of the solutions can be obtained as 
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where each diagonal element is minus the sum of the other elements in its 
column, and the following relations hold between the rates. 

{\-i)b = i[c + f + {\-Od], 

(1-0(5 + ?) = / + (i-eK 

iu^il- Or. (54) 

It is clear that for < ^ < 1, one can always find rates satisfying H54f) . With 
these rates, one has 

M = -(1+0(5 + 9), 

V ^ 9 + q^ 

e^ag + q)- (55) 



6 Concluding remarks 

Among the aims of investigating reaction-diffusion systems is to find as many n- 
point functions (of number operators) as possible. There are systems for which 
one can find certain classes of these correlators. 

The empty interval method was first introduced to investigate the probabil- 
ity of finding n neighboring empty sites, for systems consisting of one species, 
with nearest-neighbor interactions. One can generalize this method, in several 
aspects. One way is to consider systems with more-than-two-site interactions. 
Another way is to consider multi-species systems, and ask for the expectation 
of the product of certain linear combinations of the number operators. While 
this does not (generally) give the densities of each specific kind of particles, it 
does give a certain combination of these densities. 

What was introduced here, was a set of constraints the reaction-diffusion 
systems should satisfy, in order that the system be solvable through the gener- 
alization of the empty interval method. For example, in a system consisting of 
particles of one kind (with nearest-neighbor interaction), there are 12 indepen- 
dent reaction-rates. If one demands that the system be left-right symmetric, 
the number of independent rates is reduced to 7. Among these, there exists a 
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5-parameter family, the models of which are solvable through the generalized 
empty interval method. The classification of this family, and some examples, 
were discussed in section 3. An interesting problem may be to classify the 
models satisfying the solvability conditions for multi-species systems. 



7 Appendix 

We want to solve the inequalities for a given value of ^. It is seen that 

changing the signs of ^, v, and simultaneously, while keeping the signs of A, 
/X, and A fixed, does not change the inequalities. So it is sufficient to solve the 
inequalities for nonnegative ^. 

First consider the case /i = 0. If /i = 0, and ^ = 1, then one arrives at 

i) 

C = 1, ^ = 6* = 0, v<Q, 4v<X<2v. (56) 

Another case is = 0, ^ ^ 1. One can consider the subcases ^ = and ^ 7^ 0, 
and show that in both subcases, A = 9 = v = 0, which means the Hamiltonian 
is zero. So, if ^ 7^ 1, for any nontrivial solution, n 0- 

If /i 7^ 0, then it should be negative. One can divide the Hamiltonian by — /i 
(which is like taking fi = —1). The inequalities p8(l are then rewritten as 

(l-e)[-(l-0A-+2(l + r)] >0, 
(l-0[(l+e)A-+2T] >0, 
(l+e)[-(l+0A+ + 2(l-T)] >0, 
(1+0[(1-0A--2t] >0, 
(l+0'A-2e(l-r) >0, 
(l-0'A + 2e(l + T) >0, 

-(1 -e^)A + 2er + 2 > 0, (57) 

where r and A^ are defined through (|30|l . Now, two general cases occur. Either 
< ^ < 1, or 1 < ^. First, take < ^ < 1. Then, the inequalities H57|) become 

2t > -2 + (1 -e)A- =: Fi, 

2r > -(l+e)A- ^■.F2, 
:=2-(l+0A+ >2t, 
F4 (1 - 0A+ > 2r, 

2r>2-ii±il!A^:K, 



2r > -2 - -^A F^ 



2 

2T>-- + —^A^:Fr. (58) 
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To have solution for r, F3 and F4 must be greater than or equal to Fi, F2, F5, 
Fq, and F7. So, we have ten inequalities for A~ and A+. The first four, coming 
from the (F3, F4) > (Fi, F2), are 

4-(l + 0A+-(l-e)A- >0, 
2 + (l-O(A+-A-)>0, 
2-(l+O(A+-A-)>0, 

(1-0A+ + (1 + 0A" >0. (59) 

The solution to these is the interior of a tetragon ABCD. The coordinates of 
its vertices are 

^(-ItI'O. K'.^)- -('-1^ 

where the first coordinate is A~, and the second is A+. 

There remains six other inequahties for A~ and A"*", to be satisfied. As aU of 
the inequahties are hnear, it is sufficient to check the inequahties on the vertices 
of the above tetragon. In fact, we have to compare the values of F3 and F4, with 
those of -F5, Fq, and F-j, at the points A, B, C, and D. Doing so, it is seen that 
the problem is only with F^ at the point D; all other inequalities are satisfied. 
At the segments DA and CD, F4 < F3. So we have to solve the inequality 
F4 > F5. The line F4 = F5 passes through A, and intersects the segment CD 
at the point E: 

1(1-0(1 + 30' (1-0(1 + 30^' 

So, 
ii) 



< ^ < 1, ^ < (^ = -1) (A ,A+) inside the tetragon ABCF, t satisfies (|^ . 

(60) 

For the case 1 < one can still use (|57|l . But as (1 — is negative, the first 
two inequalities in l|58|) are reversed. So we have 

(Fi, F2, F3, F4) > 2r > (F5, Fg, Fj). (61) 

From these, one should have 

(Fi, F2, F3, F4) > (F5, Fg, F7), (62) 

which are twelve inequalities for A* . From Fi > Fg and F3 > F5 , one obtains 

($-l)A+-(e + l)A- >0, 
-(e-l)A+ + (e + l)A- >0, (63) 

respectively. So, one has 

(e-l)A+ = (e + l)A-=:x. (64) 
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This also means that 

Fi = = F5 = Fg = 2r. (65) 
From Fi = F^, for example, x is obtained: 

£2 - 1 

X=-^- (66) 

It is easily seen that this satisfies (|61f) . So, one arrives at 
iii) 

1<£, ^^<0{^^^-l) r^-^^, A^^^. (67) 

It is seen that here there is no allowed region for the rates, but a single point 
(apart from scaling) for each value of £. 

Finally, there remains two other special cases. First, the case £ = 0. In this 
case, the first four inequalities in H57|) are not changed, and hence the first four 
inequalities in (One should of course put £ = in them.) The vertices of 

the tetragon ABCD are now 

Ao(-l,l), Bo(l,3), Co(3,l), Do{l,-l). 

In the remaining three inequalities of (|57|l . there is no r, and one reads 

< A < 2. (68) 

It is easy to see that these are satisfied inside the square Aoi?oCo-Do. This 
square is in fact the same tetragon ABCE at the limit 1^ 0. (Note that in 
this limit, the points D and E tend to each other.) So, 
iv) 

£ = 0, 7^ (/i = — 1), (A~, A^) is inside the square Aoi?oC'o-Do- 

T satisfies (69) 

Finally, in the case £ = 1, the first two inequalities in H57|l become identities, 
and the sixth and seventh become identical to each other. So one has four 
independent inequalities: 

-A+ + 1 - r > 0, 
-r > 0, 
2A- (1 -r) > 0, 

1 + r > 0. (70) 

These give the allowed region for (A^,A+), as the region limited by the lines 
AiB' , the horizontal line passing through B' , and the line passing through Ai 
with the slope —1, containing the point (1,1), where 

Ai(0,l), S'(0,2). 
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Let us call this region S. So, 
v) 

^=1, ^^0(^ = -l), (A-,A+) is inside §. t satisfies GO)). (71) 
The cases i) to v) summarize the desired classification. 
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